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摘 要 :考虑 到 混凝土 是 一 种 非 均 质 材料 ,为 研究 新 老 混凝土 界面 籍 结 特性 和 破坏 机 理 , 利 用 
的 和 二 民生 土 黏 结 试 块 双 面 剪 切 性 能 和 增 大 截面 法 加 国 钢筋 混凝土 构件 轴 
所 六 二 和 i 人 二 入 训 红 。 结果 表明 ， 对 于 双 面 剪 切 试验 ,新 老 混 凝 土 剪 切 - 位 移 数值 曲线 与 试验 
曲线 基本 一 致 ,模型 接触 断裂 主要 发 生 在 新 老 混凝土 界面 处 ,裂缝 数量 随 着 模型 达到 峰值 后 急剧 上 
升 ,在 界面 处 形成 贯通 的 断裂 面 ;在 轴 压 试验 模拟 中 , 轴 压 力 、 混 凝 土 应 变 关系 数值 曲线 与 试验 结果 
吻合 较 好 ,数值 结果 能 较 好 地 预测 软化 阶段 ,模型 中 接触 断裂 主要 发 生 在 柱 中 部 及 下 部 位 置 , 与 实 
际 破坏 模式 相符 。 
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Abstract :Considering that concrete is a kind of heterogeneous material ,in order to study the interfacial 
bonding characteristics and failure mechanism of new and old concrete ,PFC3D mesoscopic discrete ele- 
ment software was used to simulate the double-sided shear performance of new and old concrete specimens 
and the axial compression mechanical properties of reinforced concrete members strengthened by increas- 
ing section method. The results showed that in the double-sided shear test ,the numerical curve of shear 
displacement of new and old concrete is basically consistent with the test curve ,the contact fracture of the 
model mainly occurs at the interface of new and old concrete ,and the number of cracks increases sharply 
after the model reaches the peak value ,and a transfixed fracture surface is formed at the interface. In the 
simulation of axial compression test ,the numerical curve of the relationship between axial compression and 


concrete strain is consistent with the experimental curve ,and the numerical simulation can well predict the 
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softening stage. The contact fracture mainly occurs in the middle and lower parts of the column, which is 


consistent with the actual failure mode. 


Key words :discrete element method;new and old concrete ;interface failure ;mesoscopic ; reinforcement 


关于 混凝土 的 数值 模拟 研究 已 展开 许多 研究 工 
作 , 采 用 基于 连续 介质 力学 的 有 限 元 方法 虽然 能 够 
在 一 些 条 件 下 满足 精度 要 求 ,但 不 能 描述 颗粒 间 的 
复杂 相互 作用 及 高 度 非 线 性 行为 , 亦 无 法 正确 刻画 
材料 的 流动 变形 特征 。 

颗粒 流 (particle flow code ,PFC ) 031 是 一 种 细 观 
离散 元 软件 , 它 通过 模拟 圆 形 颗粒 介质 的 运动 及 相 
互 作用 来 研究 颗粒 介质 特性 。PFC 通过 颗粒 单元 构 
造 试 样 几何 形状 .接触 本 构 关 系 形成 相互 作用 、 和 迭代 
分 析 法 更 新 颗粒 的 状态 ,最 终 使 得 数值 试 样 的 宏观 
力学 特性 逼近 真实 材料 的 力学 特性 ,并 且 呈 现 可 视 
修 的 结果 。PFC 基于 非 连续 介质 力学 理论 设计 ,有 
着 道 合 混凝土 材料 的 本 构 接触 模型 ,并 且 对 介质 的 
:小 分 离 等 非 连续 现象 有 较 好 的 模拟 效果 ,在 混 凝 
长 售 域 有 着 较 好 的 应 用 前 景 。 张 正 至 等 5 模拟 了 二 
维 混 凝 土 单 铀 压缩 破坏 全 过 程 , 探讨 了 数值 试验 中 
应 注 -应 变 曲线 及 体积 应 变 的 变化 规律 ,分 析 了 颗粒 
怀 请 颗粒 数目 等 细 观 参数 对 模拟 结果 的 影响 ; 许 尚 
术 等 “进行 了 单 级 配 和 全 级 配 单 轴 压 缩 试验 二 维 离 
散 过 模拟 ,得 出 骨 料 空间 位 置 和 数量 ,大 小 均 会 影响 
混凝土 的 抗 压 性 能 ; 赵 全 满 等 ”利用 PFC2D 对 混 北 
隆 搞 压强 度 进行 了 细 观 数值 模拟 ,分 析 细 观 参数 对 
弹性 模 量 和 抗 压 强度 的 影响 ,通过 校 核 确定 细 观 参 


炜 全 上 9 利用 PFC2D 探究 双 轴 压缩 过 程 中 颗粒 物质 
接触 力 和 力 链 特性 ,通过 接触 力 大 小 、 角 度 分 布 规律 
等 细 观 层次 分 析 来 实现 对 力 链 方向 性 的 量化 描述 。 
上 述 研 究 均 基 于 二 维 离散 元 模拟 ;但 实际 混凝土 的 
裂 颖 发生、 扩展 都 是 三 维 的 ,车 将 混凝土 试 件 破 坏 过 
程 看 成 是 二 维 平面 问题 , 则 不 能 真实 反映 混凝土 试 
件 复杂 的 应 力 应 变 关系 。 此 外 ,通过 有 限 元 方法 解 
决 新 老 混凝土 界面 许 结 "问题 并 不 理想 ,例如 混 凝 
土 与 混凝土 界面 在 剪 力 和 弯 矩 组 合作 用 下 的 力学 性 
能 和 破坏 形态 .新 旧 混 凝 土 错 固 钢筋 混凝土 模型 的 
破坏 模式 、 界 面 剪 切 滑 移 等 ,而 细 观 离散 元 方法 为 该 
问题 提供 了 解决 的 可 能 性 ;但 是 目前 相关 的 研究 依 
然 匮乏 。 

文献 [10-11] 分 别 利 用 PFC3D 软件 研究 了 不 同 
孔隙 率 下 无 砂 混凝土 单 轴 压 缩 下 的 受 力 及 破坏 情况 


和 端 部 摩擦 对 不 同形 状 混凝土 抗 压强 度 、 变 形 和 破 
坏 模式 的 影响 。 本 研究 首先 采用 PFC3D 离散 元 软 
件 对 新 老 混凝土 笑 结 试 块 双 面 前 切 试验 进行 数值 模 
拟 ,并 与 试验 数据 进行 对 比 ,从 细 观 角度 研究 其 力学 
特性 裂缝 演化 过 程 和 分 布 规律 , 以 验证 离散 元 方法 
对 于 界面 问题 的 适用 性 。 然 后 ,对 增 大 截面 法 加 固 
钢筋 混凝土 轴 压 试验 进行 数值 模拟 ,并 与 现 有 数据 
进行 对 比 ,分 析 了 新 老 混 凝 土 界面 的 受 力 特性 与 破 
坏 模式 。 本 研究 研究 成 果 将 对 离散 元 方法 在 混凝土 
结构 的 加 固 应 用 提供 一 定 的 参考 。 


1 PFC3D 离散 元 数值 模拟 


接触 本 构 是 离散 元 计算 方法 的 关键 ,能 反映 微 
细 观 介质 的 力学 特性 。 其 中 平行 茜 结 模型 适用 于 系 
结 材 料 的 力学 行为 (类似 于 环 氧 树脂 胶合 玻璃 珠 、 水 
泥 间 上 骨 料 的 灰 结 ) ” ,可 视 作 一 组 弹簧 ,均匀 分 布 于 
接触 面 与 接触 点 。 这 些 弹簧 与 线性 元 件 弹 簧 平行 ， 
平行 键 产生 后 ,接触 处 发 生 相 对 运动 ,颗粒 具有 滑动 
的 可 能 性 ,平行 务 结 可 以 在 颗粒 之 间 传 递 力 和 力 
和 矩 ” ; 当 法 向 应 力 超过 其 相应 的 务 结 强度 ,平行 秋 
结 断 裂 , 模 型 退化 为 线性 模型 ,此 时 相互 接触 颗粒 蔚 
结 发 生 断 裂 形 成 微 损伤 ,进而 扩展 成 破裂 面 使 试 样 
破坏 ,这 与 实际 混凝土 材料 的 损伤 发 展 相符 ” 。 平 
行 务 结 接触 如 图 1 所 示 。 
抗 转动 生 结 


抗 前 切 条 结 猎 结 刚度 


番 结 刚度 破坏 
图 1 平行 务 结 接触 示意 图 
Fig.1 Parallel bond contact 
平行 黏 结 中 的 力 -位 移 法 则 满足 
F.=F+F +F (1) 


M, = 站 (2) 
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式 中 :Fr 为 线 弹 性 力 ; Fr" 为 阻尼 力 ;为 平行 黏 结 
力 ; 政 为 平行 竺 结 弯 矩 。 其 中 , 线 弹性 力 和 阻尼 力 在 
线性 模型 中 更 新 ,而 平行 香 结 力 和 弯 矩 在 平行 香 结 
模型 中 更 新 。 
平行 竺 结 力 分 解 为 法 向 力 和 切 向 力 ,而 平行 符 
结 弯 矩 由 扭矩 和 弯 矩 合成 。 
F= -Fn 二 下 (3) 
WH=Mh, +W, (4) 
式 中 ,>0 为 拉力 。 平 行 黏 结 剪 力 和 弯 算 在 坐标 系 
可 分 解 为 
F = 下 +h,i, (5) 
M, = M,s, + M,i. (6) 
当 平 行 茜 结 的 联系 通过 黏 结 方法 被 建立 ,平行 
符 结 力 和 过 和 矩 初始 值 为 0。 平 行 黏 结 模型 在 两 颗 球 
体 表 面 提供 相互 作用 , 当 黏 结 发 生 破坏 时 ,相互 作用 
Ce 
例 移 增 量 积累 。 
已 5 = A, (7) 
或 册 A5, 为 法 向 位 移 增 量 。 
© 


2 中 双 面 剪 切 试验 


©O 
CN 双 面 剪 切 试验 采用 2/3 体积 老 混凝土 与 1/3 体 
积 新 混凝土 双 面 黏 结 的 边 长 100 mm 复合 立方 体 斌 
块 x 处 验 试 块 加 载 情 况 如 图 2 所 示 。 老 混凝土 采用 龄 
期 GR0 d 边 长 100 mm 的 C30 立方 体 试 块 。 将 试 块 平均 
分 段 3 等 份 。 对 切割 后 的 老 混凝土 试 块 采取 人 工 羡 毛 
处 理 , 玄 去 备 结 面 的 浮 桨 和 部 分 水 泥 , 形 成 四 凸 不 平 表 
面 ;如 图 3 所 示 。 新 增 混凝土 强度 等 级 为 C35, 共 制备 
16 组 48 个 试 块 ,标准 养护 28 d。 试验 在 WE-300 型 万 
能 试验 机 上 进行 ,加 载 速度 为 0.05 MPavsi5 1 。 
;人 
| 4 医 结 面 
3 新 混凝土 
1.2 老 混凝土 


Fig.2 Double-sided shear test 
利用 Fish 语言 和 PFC 内 置 命令 制作 剪 切 盒 ， 
ball 颗粒 填充 后 删除 新 老 混凝土 交界 wall ,产生 接触 
条 件 后 静 置 收敛 ;随后 通过 伺服 控制 促使 颗粒 位 置 
不 断 调整 ,实现 区 域 均匀 孔 际 率 、 均 匀 应 力 ; 再 赋予 
接触 参数 ,加 载 方式 以 速度 加 载 控 制 , 数 值 模型 如 图 
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4 所 示 。 其 中 新 老 混凝土 界面 粗糙 度 处 理 如 图 5 所 
示 。 建 模 时 首先 对 混凝土 进行 分 组 , wall 左 侧 为 老 
混凝土 , 右 侧 为 新 混凝土 。 在 交界 wall 移 除 之 前 , 颗 
粒 有 一 定 的 重 辣 量 , 紧 贴 合 于 wall。 移 除 wall, 此 时 
颗粒 会 弹 开 ,计算 至 平衡 ,使 得 原本 界面 处 位 置 填充 
满 颗 粒 。 此 时 , 原 wall 处 新 老 混凝土 颗粒 之 间 有 了 
接触 条 件 ,界面 分 布 不 规则 曲线 。 新 老 混凝土 颗粒 
之 间 颗 粒 采 用 线性 接触 本 构 。 对 于 ball 单元 生成 的 
模型 ,通过 建 模 前 预 设 控制 参数 emod 来 控制 粗糙 度 
的 大 小 ,数值 越 大 黏 结 越 紧 密 , 说 明 越 粗糙 ,本 研究 
采用 1 GPa。 


到 3 老 混凝土 凿 毛 处 理 


Fig.3 Roughening of old concrete 


图 4 双 面 剪 切 数值 模拟 示意 图 


Fig.4 Double-sided shear numerical simulation 


墙 体 删除 墙 体 


图 5 界面 粗糙 度 示 意图 


Fig.5 Interface roughness 
2.1 细 观 力学 参数 的 选取 
PFC 的 基本 思想 在 很 大 程度 上 不 同 于 其 他 连续 


和 非 连续 力学 理论 方法 ,导致 PFC 与 常规 数值 方法 
之 间 存 在 一 定 的 差异 ,模型 介质 的 宏观 基本 物理 力 
学 参数 不 能 直接 与 颗粒 细 观 的 结构 参数 联系 ,需要 
通过 数值 试验 标定 程序 得 到 与 物理 模型 宏观 力学 相 
对 应 的 细 观 参数 。 不 同 于 以 往 直 接 对 接触 刚度 和 
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上 赋值 ,由 于 颗粒 形状 和 粒 径 的 差异 ,采用 有 效 模 量 
emod 对 不 同 接触 的 黏 结 自动 计算 其 刚度 。 利 用 平 
行 黏 结 线性 快速 标定 法 ,得 出 的 模型 细 观 参数 如 表 
1 所 示 。 


表 1 双 面 剪 切 试 件 细 观 参数 


Tab.1 Meso parameters of double-sided shear test 


混凝土 ”有 效 模 量 /GPa ”刚度 比 。 平行 笑 结 有 效 模 量 /GPa 


平行 儿 结 法 向 强度 /MPa ”平行 香 结 切 向 强度 /MPa ”平行 寿 结 刚度 比 


老 混凝土 5.4 1 5.4 
新 混凝土 2.5 1 2.5 
2.2 前 切 滑 移 曲线 


图 6 为 新 老 混凝土 竺 结 试 块 双 面 剪 切 试验 曲线 
与 数值 曲线 的 比较 。 从 图 6 中 可 以 看 出 ,试验 曲线 
了 5 移 为 0.25 mm 之 前 呈 线 性 增长 趋势 ,之 后 曲线 
斜 痢 变 小 ,强度 逐渐 上 升 ,达到 峰值 后 进入 软化 阶 
段 3 而 模拟 曲线 在 达到 峰值 前 基本 呈 线 性 增长 ,达到 
后 开始 下 降 , 相 较 于 试验 曲线 软化 阶段 更 为 明 
最 * 昌 然 在 位 移 到 0.3 mm 后 ,模拟 剪 切 强度 与 试验 
尊 园 强度 具有 一 定 差异 ,但 是 从 整体 上 看 ,两 者 差异 
,模拟 效果 较为 理想 。 


20 视 


一 一 模拟 
一 一 试验 
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位 移 /mm 


图 6 剪 切 滑 移 对 比 
Fig.6 Comparison of Shear-slip 


图 7 为 双 面 剪 切 试 块 的 位 移 云图 。 从 图 7 中 可 


38 38 1 


42 40 1 


部 深蓝 色 区 域 时 ,位 移 量 差 值 增 大 ,说 明 此 时 裂缝 出 
现 的 位 置 会 围绕 着 右 侧 深蓝 色 区 域 。 


颗粒 位 移 
颗粒 (20 663) 
9.049 1X10” 
9.0000X10” 
8.500 OX 107 
8.000 0X107 
7.5000X10” 
7.0000X10” 
6.500 0X10” 
6.0000X10? 
5.5000X10” 
5.000 0X10” 
4.5000X10” 
4.0000X10” 
3.5000X10” 
3.0000X10” 
2.500 0X10” 
2.0000X10” 
1.5000X10” 
1.0000X10” 
5.0000X10” 
1.555 8X10° 


图 7 双 前 试 块 位 移 云 图 
Fig.7 Displacement nephogram of 
double-sided shear test block 


2.3 ”损伤 断裂 演化 


双 面 剪 切 模型 裂隙 演化 过 程 如 图 8 所 示 。 当 时 
步 为 5 000 时 ,新 老 混凝土 上 表面 交界 处 及 试 块 四 周 
出 现 少 量 的 红色 圆 片 状 裂 际 网 格 ,数量 为 341; 当时 
步 为 10 000 时 , 裂 际 网 格 集中 于 加 载 界面 处 ,数量 为 
1181 ,说 明 加 载 处 界面 优先 出 现 破 坏 裂 颖 ; 当时 步 
为 20 000 时 ,整个 试 块 中 部 以 上 的 界面 处 已 布 满 密 


以 看 出 ,新 混凝土 上 1/3 部 分 毗邻 加 载 面 ,位 移 量 最 
大 ;中 下 部 分 位 移 量 递减 ,底部 左右 两 侧 深 蓝 色 区 域 
位 移 量 最 小 。 两 侧 老 混凝土 整体 位 移 量 差 较 小 ,位 
移 较 大 人 处 集中 于 中 部 ;界面 处 由 于 新 老 混凝土 的 位 
移 量 差 值 的 大 小 而 表现 出 不 同 的 颗粒 颜色 分 布 , 界 
面 上 部 相差 最 大 ,也 是 初始 裂缝 最 容易 出 现 的 地 方 。 
位 移 量 差 值 自 上 而 下 逐渐 减 小 ,直至 下 1/3 段 时 , 右 
侧 界 面 位 移 量 差 值 几 乎 为 0, 和 裂缝 不易 产生 ;接近 底 


集 的 裂 孙 网 格 , 数 量 为 4490, 此 时 裂缝 已 进一步 扩 
大 ;当时 步 为 30000 时 , 左 侧 交 界面 处 已 完全 布 满 裂 
隙 网 格 ,而 右 侧 交界 面 上 2/3 段 裂缝 沿 着 竖 直 方向 
拓展 ,下 1/3 段 裂缝 向 老 混凝土 一 侧 偏 移 , 最 后 又 折 
返 于 新 老 混凝土 交界 处 底部 , 裂 际 数量 为 12 309。 
从 整体 上 看 ,裂缝 的 形成 .扩展 发 生 在 新 老 混凝土 界 
面 处 (1 ,说 明 PFC3D 对 新 老 混凝土 双 面 剪 切 界 面 
问题 的 研究 是 可 行 的 。 
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QJ 图 9 为 试 件 双 面 剪 切 破坏 图 。 从 图 9 中 可 以 看 
通 裂 颖 主要 发 生 于 新 老 混凝土 界面 处 , 试 块 表 
恒 深 有 少许 的 混凝土 脱落 。 当 试验 机 加 载 至 峰值 荷 
,伴随 着 一 声 脆 响 , 秋 结 面 出 现 裂 锋 , 试 件 随机 
分 离 。 破 坏 面 较为 平坦, 属于 与 型 的 脆性 破坏 ， 
相思 于 数值 模拟 ,破坏 形态 大 致 吻合 。 


图 9 双 面 剪 切 破坏 示意 图 
Fig.9 Diagram of double side shear failure 


3 ” 增 大 截面 法 加 固 钢筋 混凝土 构件 轴 
压 数值 模拟 


为 进一步 研究 不 同 加 载 方式 下 PFC3D 对 于 界 
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(d) 时 步 30 000 


图 8 双 面 剪 切 模型 裂隙 演 化 过 程 


Fig.8 Fracture damage evolution of double-sided shear model 


面 研究 方向 的 可 行 性 ,对 增 大 截面 法 加 固 钢筋 混 凝 
土 轴 压 构件 进行 数值 模拟 。 试 验 原 柱 截面 "” 尺 十 
200 mm x200 mm , 高 1 000 mm; 加 固 后 截面 尺寸 扩 
大 到 320 mm x 320 mm; 纵 筋 采用 直径 14 mm 的 
HRB335 级 钢筋 , 短 筋 采用 直径 6 mm 的 HPB300 级 
钢筋 ,间距 150 mm; 加 固 前 首先 对 原 柱 进行 溺 毛 处 
理 , 新 老 混凝土 界面 用 水 泥 净 浆 进行 处 理 。 数 值 模 
拟 中 建立 内 外 墙 体 ;采用 Clump 命令 填充 ,与 传统 的 
球体 颗粒 填充 相 比 ,通过 外 部 导入 方式 引进 多 种 几 
何 形状 (Ceometry 命令 ) 以 模拟 不 规则 骨 料 ,更 具 真 
实 性 ,如 图 10 所 示 ; 利 用 Fish 语言 命令 来 实现 纵 筋 
和 短 筋 的 建 模 ; 顶 面 和 底面 为 刚性 墙 体 以 模拟 加 载 
板 , 内 外 侧面 均 为 自由 边界 。 钢 筋 混 北 土 整 柱 数值 
模型 如 图 11 所 示 。 


602 人 9 
© 


图 10 骨 料 仿真 样式 


Fig. 10 Style of aggregate simulation 


混凝土 
| 四 骨 料 
| 国 颗粒 半径 
| 1.200 0X103 
国 11500X103 
| 二 
dd 。 m 1.0000X102 
9.500 0% 103 
D | 90000X103 
1 诗 = 认 有 
a 2 
/由 = 由 i 
人 [| 5 500 0X 103 
p 让 
Flt 4.000 0X 103 
CC 让 
sr 2 500 0X 103 
233797X10° 
钢筋 的 布置 颗粒 尺寸 
(a) 整体 组 成 
™ 
国 
©O (b) 局 部 组 成 
CA 图 11， 增 大 鹤 面 法 加 固 钢筋 混凝土 柱 数值 模型 
Cig.11 Numerical model of reinforced concrete column 
GAN strengthened by enlarging section method 


三 与 上 述 双 剪 数值 试验 建立 单元 不 同 ,采用 
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Clump 单元 填充 。 因 其 不 规则 的 外 形 ,所 以 界面 在 
模型 加 载 前 的 初始 收敛 阶段 后 会 有 一 定 的 四 凸 ,类 
似 于 实际 中 的 浅 毛 处 理 。 界 面 处 进行 界面 接触 分 
区 ,接触 本 构 采 用 平行 黏 结 模型 ,填充 新 老 混 凝 土 界 
面 粗糙 度 处 理 如 图 12 所 示 。 为 了 实现 对 模型 中 粗 
糙 度 的 控制 ,通过 改变 模板 填充 Bubplepack 命令 中 
的 Distance(0 ~ 180) 和 Ratio(0 ~ 1) 。 考 虑 计算 效 
率 和 硬件 设施 条 件 , 这 里 将 其 分 别 设 为 120 和 0.3。 


删除 墙 体 
0 


图 12 界面 粗糙 度 示意 图 


Fig. 12 Interface roughness 
3.1 细 观 力学 参数 选取 
由 于 接触 分 组 较 多 ,赋予 接触 参数 较 多 ,参数 之 
间 也 相互 影响 ,平行 务 结 快速 标定 法 不 能 完全 解决 ， 
参考 Fakhimi 等 "提出 的 量 纲 分 析 法 ,结合 快速 标 
定 法 ,最 后 得 出 模型 细 观 参数 , 见 表 2。 


>< 表 2 钢筋 混凝土 细 观 参数 
© 


Tab.2 Meso parameters of reinforced concrete 


有 效 模 量 /GPa 


刚度 比 。 平行 莫 结 有 效 模 量 /GPa 


平行 灰 结 法 向 强度 AMPa ”平行 香 结 切 向 强度 A/MPa ”平行 寿 结 刚度 比 


wk 1 1 


旧 柱 骨 料 5 1 10 
新 柱 砂 浆 1.7 1 2 
新 柱 骨 料 5 1 10 
钢筋 10 1 10 


3.2 轴 压 力 与 应 变 曲线 


对 混凝土 柱 上 下 刚性 加 载 板 使 用 速度 加 载 法 ， 
利用 Fish 语言 对 轴 压 力 和 混凝土 应 变 关系 进行 数 
据 提 取 ,最 终 将 模拟 结果 与 试验 结果 '” 进行 对 比 ， 
结果 如 图 13 所 示 。 从 图 13 中 可 以 看 出 ,在 曲线 的 
上 升 阶段 ,新 老 混凝土 数值 模拟 结果 与 试验 结果 基 
本 一 致 ,此 时 数值 模拟 中 新 老 混 凝 土 曲线 也 比较 吻 
合 ,说 明 两 者 在 加 载 初 期 能 较 好 地 协同 工作 。 在 下 


40 40 1 
200 200 1 
118 118 1 
200 200 1 
500 500 1 


降 阶 段 ,模拟 段 曲 线 呈 凸 型 ,塑性 应 变 不 断 增加 , 软 
化 现象 较为 明显 ;在 模拟 曲线 上 升 中 后 段 和 下 降 段 ， 
新 老 混凝土 模拟 曲线 有 略微 差异 ,说 明 新 老 混凝土 
在 相同 轴 压 力 下 会 有 一 定 的 相对 位 移 。 这 是 因为 新 
老 混凝土 黏 结 面 为 薄弱 面 。 


3.3 ”损伤 断裂 演化 


混凝土 损伤 裂隙 演化 过 程 如 图 14 所 示 。 通 过 
PFC 中 的 内 置 Fish 语言 编程 来 实现 对 颗粒 间接 触 断 
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裂 的 实时 监测 ,其 中 DFN (discrete fracture network ) 
为 离散 裂 际 网 格 , 当 颗 粒 之 间 的 接触 力 大 于 黏 结 强 
度 时 ,就 会 生成 圆 盘 状 ( 二 维 为 直线 ) ,说 明 黏 结 材 
料 接触 断裂 。 当 时 步 为 4000 时 , 裂 际 均 匀 分 布 于 混 
凝 土 四 周 ; 随 着 时 步 的 增加 , 裂 际 逐 渐 密集 ,此 时 和 裂 
际 直 径 较 小 ,损伤 较 小 ,如 图 14(b) 所 示 ; 当 运 行 至 
7 000 步 时 ,混凝土 中 部 、 上 部 和 下 部 均 有 较 大 裂 院 
出 现 ,但 是 主要 集中 于 中 下 部 ;当时 步 为 8000 时 , 混 
凝 土 中 部 内 外 侧面 均 出 现 较 大 和 裂 际 ,下 部 裂 际 主 要 
集中 于 内 部 ; 当时 步 为 9 000 时 ,混凝土 中 部 外 表面 
存在 明显 的 主 裂 际 ,说 明 外 表面 有 明显 的 混凝土 脱 
落 , 下 部 裂 际 也 愈 发 密集 ,破坏 较为 明显 。 
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图 13 混凝土 轴 压 力 与 应 变 曲线 


Fig. 13 ”Axial compression and strain curve of concrete 
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图 14 


混凝土 数值 模型 裂 孙 演 化 过 程 


Fig. 14 Fracture damage evolution of numerical concrete model 


4 结 论 


1) 对 于 双 面 剪 切 试验 ,模拟 得 到 的 剪 切 滑 移 数 
值 曲 线 与 试验 曲线 能 较 好 地 吻合 ; 随 着 加 载 力 的 逐 
渐 增 大 ,初始 裂 颖 出 现在 了 加 载 面 处 的 新 老 混 凝 土 
交界 位 置 , 裂 颖 沿 着 交界 面 逐 渐 延 伸 ; 在 模型 达到 峰 
值 强 度 后 ， 人 康 率 急剧 增加 ;最 终 裂 颖 布 满 交 
界面 ,形成 贯通 的 断裂 面 ,说明 新 老 混 凝 土 交界 面 是 
薄弱 面 。 


2) 利 用 Clump 命令 进行 仿真 模拟 ,保证 与 实际 
骨 料 形态 特征 一 致 , 增 大 颗粒 数量 来 提高 模拟 精度 。 
轴 压 试验 中 ,将 得 出 的 轴 压 力 ` 混 凝 土 应 变 关系 曲线 
与 观 有 试验 曲线 比较 , 抗 压 强度 一 致 ,并 且 模拟 曲线 
给 出 了 软化 阶段 ,说 明了 本 研究 方法 对 混凝土 力学 
特性 预测 的 可 行 性 。 研 究 发 现 ,裂缝 集中 于 柱 中 下 
部 ,中 部 内 外 表面 均 有 较 大 裂缝 , 受 损 严重 ; 柱 下 部 
裂 际 较 小 ,但 密集 分 布 于 内 部 ,与 试验 破坏 模式 较为 


Ss: 
Bs 


3) PFC 软件 对 混凝土 力学 性 能 \ 破 坏 形态 


和 弄 


缝 拓展 过 程 的 追踪 有 着 良好 的 模拟 效果 ,能 有 效 地 
预测 加 固 混凝土 起 始 裂 颖 的 位 置 及 主 裂 缝 的 形态 ， 
因此 可 以 针对 起 裂 点 处 进行 二 次 加 固 工作 。 
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